Geometric fluctuations in a two-dimensional quantum antiferromagnet 
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We consider the effects of random fluctuations in the local geometry on the ground state properties 
of a two-dimensional quantum antiferromagnet. We analyse the behavior of spins described by the 
Heisenberg model as a function of what we call phason disorder, following a terminology used for 
aperiodic systems. The calculations were carried out both within linear spin wave theory and using 
quantum Monte Carlo simulations. An "order by disorder" phenomenon is observed in this model, 
wherein antiferromagnetism is found to be enhanced by phason disorder. The value of the staggered 
order parameter increases with the number of defects, accompanied by an increase in the ground 
state energy of the system. 
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Antiferromagnetic Heisenberg models are particularly 
interesting in the case of two dimensions, where quan- 
tum corrections have been studied in a number of lat- 
tices, with and without disorder 0, Q. For clean sys- 
tems, and in the absence of frustrating interactions, the 
ground state is believed to possess collinear Neel semi- 
classical order even in the extreme quantum limit of spin 
S = 1/2. Disordered two dimensional Heisenberg mod- 
els have also been much studied in the literature. These 
include models of random exchange, in which the spin ex- 
change couplings are taken from some distribution 0, Q 
and site dilution models in which a fraction of sites are 
unoccupied by spins 0-0] ■ I n tms paper, we consider 
disorder of geometrical origin, in which random fluctua- 
tions cause a local switching of the connections between 
the sites. In quasicrystals, this type of fluctuation of the 
local environments is well known by the name of "phason 
flip disorder" , and has been proposed as one of the im- 
portant mechanisms for the formation of such structures 
based on their entropy § . Such defects could in principle 
exist in regular crystals as well, and in this paper we ad- 
dress the question of their consequences for the quantum 
antiferromagnetic Heisenberg model. This type of dis- 
order does not affect the classical energy of the system 
because the overall number of bonds and the strength of 
the couplings is not affected by the randomness. Quan- 
tum fluctuations are modified, however, and we will see 
that a type of order-by-disorder phenomenon Q occurs 
in this model. 

In contrast to the exchange- and site-disorder models, 
where antiferromagnetism is weakened and ultimately 
destroyed for sufficiently strong disorder, phason disor- 
der strengthens the order parameter in our model. The 
ground state energy is increased by defects, and the in- 
teraction between defects. 

The lattice we consider is the staggered dice lattice, 




FIG. 1: The staggered dice lattice without disorder (left) 
and with phason disorder (right). 



closely related to the usual dice, or T3, lattice. The lat- 
ter has been investigated for unusual electronic proper- 
ties associated with its flat band and Dirac cones [Ill - ll3l | . 
These lattices correspond to different decorations of the 
honeycomb lattice, as we will now explain. The staggered 
dice lattice without disorder is shown in the left hand 
panel of FigJT] The bonds corresponding to the honey- 
comb lattice are shown with thick lines, while the bonds 
in the "interiors" of the hexagons are shown by thin lines. 
The sites of the hexagonal lattice form hexagonal cages 
which enclose a site in the interior. The "interior" sites 
are connected to three of the cage sites in one of two pos- 
sible ways: one can draw three bonds that have the form 
of a "Y" that faces either upwards (U) or downwards 
(D). In the ordinary dice (T3) lattice, all hexagons have 
identical interiors, whereas in the staggered dice lattice, 
there is an alternation of U and D, along one of the two 
crystallographic directions. A phason flip, in this system, 
corresponds to a local fluctuation in which a U-site trans- 
forms to a D-site or vice- versa (see Fig|2]). When phason 
flips are randomly introduced in the SDL, one obtains 
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FIG. 2: Bond configurations before (left, in the U configura- 
tion) and after (right, in the D configuration) a single phason 
flip. 



a structure such as the one illustrated in the right hand 
figure of FigfTJ The pure dice lattice (i.e. without dis- 
order) can be realized in optical lattices, as discussed in 
[rij . and it would be interesting to similarly realize the 
phason disordered model experimentally. 

In the staggered dice lattice phason flip disorder can- 
not be increased to arbitrarily high values. In fact, the 
problem shows a symmetry with respect to the point 
A = i, where < A < 1 is the concentration of phason 
flips. Linear spin wave theory is used to study the sys- 
tem within the entire range of allowed values, and it gives 
results in agreement with the results obtained by quan- 
tum Monte Carlo calculation using the stochastic series 
expansion (SSE) method (T3 |. 

The dice and staggered dice lattices are both bipar- 
tite, so that nearest-neighbor Heisenberg exchange inter- 
actions are unfrustrated. However, there is a significant 
difference between these models. The dice lattice has two 
in-equivalent sublattices, with a different number of spins 
on each sublattice Na = 2Nb where A and B denote the 
two sublattices, and the ground state of the Heisenberg 
model is thus a ferrimagnet (some properties of which 
were investigated in Ref. 15). The ground state of the an- 
tifcrromagnetic quantum Heisenberg model on the stag- 
gered dice lattice of Fig. QJi, on the other hand, is an an- 
tifcrromagnet with two identical sublattices (Na = Nb)- 
The structural and magnetic unit cells are identical, con- 
sisting of six sites, two of which correspond to centers of 
hexagons. The coordination number z = 3 for the inte- 
rior sites, while the "cage" sites have z = 4 and z = 5. 
The number of sites of each type is equal to N/3, and 
the average coordination number on the staggered dice 
lattice z = 4. 

Consider now the consequences of a phason flip at one 
of the interior sites that changes the local bond configura- 
tion from U to D. In this process, the central site retains 
the coordination number z = 3 while the three previous 
neighbors, belonging to the surrounding hexagon, are re- 
placed by three new ones. As can be seen in Figj2l the 
interior site changes its sublattice after such a flip. In 
order to conserve the equality Na = Nb, we therefore 
perform the same number of such phason flips on each of 
the sublattices. When all possible sites are flipped, i.e., 



all the U-sites are changed to D-sites and vice-versa, the 
original lattice is recovered upto a global shift along the 
y-direction (cf. Fig. |TJl) . Thus, the maximum number 
N m ax of possible flips Nfn ps is equal to the number of 
interior sites, N max — N/3, where N denotes the total 
number of lattice sites. If we define a disorder strength 
as A = N flips /Nmax, the problem is symmetric with re- 
spect to A = 1/2. 

The specific model that we consider is the spin-1/2 
quantum antiferromagnetic Heisenberg model, with the 
Hamiltonian 



(1) 



where Si are the spin-1/2 operators for the sites i of the 
staggered dice lattice, and where we assume that the cou- 
plings J > are non-zero for all pairs of spins i,j which 
are joined by an edge. Since this system is unfrustrated, 
one expects that the SU(2) symmetry will be broken in 
the ground state, which will have Neel order in which 
spins are aligned along a direction, say, the z-axis. The 
local staggered magnetization m si denotes the absolute 
value of the component of the local magnetic moment 
in the z-direction (the direction of symmetry breaking), 
m sl = \(Siz)\. 

A linearized theory is obtained using the Holstein- 
Primakoff transformation [l6| and keeping only terms 
quadratic in the bosonic operators of the resulting Hamil- 
tonian. The ground state energy per site and the lo- 
cal staggered magnetizations were computed numerically 
for the disordered lattices, as well as for clean sam- 
ples 
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18j . In the latter case one can also reduce the 
Hamiltonian to a 6 by 6 matrix using the Bloch theorem, 
and solve analytically at least in the small k limit. The 
linearized Hamiltonian reads 



Hlsw = ~JS(S + l)N b + jsn 2 , 
n 2 = J2(a\ai + bjb] + a\b] + bjOi) 



(2) 



where the aj and bj are the boson destruction operators 
for sites on the A- and B-sublattice respectively. The 
number of bonds Nb = 2N. After the diagonalization, 
one obtains 



Hlsw = NE Q + ^^l\al^ 
E = -JS(S+l)N b /N + JSY J u»- (3) 



where the index a = ± labels the two families of 
modes that are obtained for this bipartite lattice, and 



H= 1,2, 



JV 



— 1 . Eq is the ground state energy per site. 



The numerical diagonalizations of the LSW Hamiltonian 
were carried out for several system sizes upto a maximum 
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size of N = 4800. For the clean staggered dice lattice, we 
find that the ground state energy per site, extrapolated 
to the thermodynamic limit, gives E = — 0.65f7(7)J 
within LSW theory, which compares well to the result 
Eq = —0.6639(1) J obtained from QMC simulations. 

We now consider the changes due to phason flips. Con- 
sider first the system with two flips, one on each sub- 
lattice. The new Hamiltonian can be written H'=H+ 
AH, where AH is the change due to a flip on the A- 
sublattice at position r\ and one on the B-sublattice at 
position r*2- The interaction energy between two phason 
defects can be calculated to second order in AH, and a 
naive power counting argument indicates that the inter- 
action should decay as 1/r 3 at large distances [19| (where 
r = |r*i —7%|). Our numerical results for the largest sam- 
ple studied (N = 4800) are consistent with a 1/r 3 decay 
at large distances, although bigger samples are needed in 
order to accurately determine the asymptotic behavior. 
The interaction, which is attractive at short range, tends 
at large distances to a value 2Efu p , where Efu p ps 0.057J 
is the energy of creating a single flip. This value can be 
compared with, for example, the energy of eliminating a 
single spin, which is about 0.58 J in LSW theory [20j. It 
is interesting to compare with results obtained in Ref. [21 
for the interaction between two magnetic defects in the 
square lattice antiferromagnet. The authors carried out 
a perturbation expansion in g, the ratio of impurity spin 
coupling J' to the pure coupling J, and found that the 
long range interaction between the defects was indepen- 
dent of the coupling and decays as 1/r 3 . The sign of the 
interaction was found to be negative for defects on oppo- 
site sublattices, and positive for defects on the same sub- 
lattice. We find similarly that the phason-flip interaction 
is attractive (repulsive) for phasons on different (same) 
sublattice. Work on the two flip problem is in progress 
and will be reported elsewhere. 

We consider now the evolution of ground state Eq(A) 
for finite values of the parameter < A < 1 . The ground 
state energy per spin as calculated by LSW is shown 
in Fig. |3k, while Fig. [3jj shows the comparison between 
LSW theory and QM in the low-A regime (system size 
N = 1200). The ground state energy rises at first lin- 
early, before the curve begins to fall over with increasing 
A due to interactions between phason defects. As ex- 
pected, the curve is symmetric with respect to A = 1/2. 

This behavior is well fitted by the expression 



E (A) = E + A(E max - E )(A - A 2 ) 



(4) 



with E Q = -0.6514J and E max = -0.6447 J. The initial 
linear dependence of Eo(A) is due to the energy of for- 
mation of phasons, whose interactions are negligible in 
the dilute limit. The subsequent curvature of the energy 
and its quadratic dependence on A can be explained by 
the fact that the interaction between defects decay very 
rapidly with distance. This leads to an energy contribu- 
tion proportional to the square of the number of phason 
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FIG. 3: a) The dependence of the ground state energy Eo(A) 
on the disorder strength A calculated within LSW theory, 
including a quadratic fit (cf. text), b) Comparison between 
LSW and QMC results for the range < A < 0.2. Here, 
N = 1200. 



flips. Since the disorder is bounded and there is symme- 
try around A = 1/2, this leads to the behavior shown in 
Fig. El The phason disorder results in diminishing quan- 
tum fluctuations - the ground state energy rises towards 
the classical value of E c i = —0.5 J. This behavior is ex- 
plained by the evolution of the density of states with in- 
creasing disorder. The integrated density of states N(E) 
for a sample of 1200 spins is shown in FigJU for the pure 
system (blue curve) and after averaging over samples for 
a fixed value A = 0.025 (red curve). The principal ef- 
fects of the disorder are to i) fill in the gaps and ii) intro- 
duce high energy states around E ~ 2.75 J (states arising 
from the new environments created by disorder, namely 
sites with z=6). The result is an increase of the integral 
over energies, which enters in the expression for the LSW 
ground state energy (Eqj3j), explaining the observed rise 
in Eq with increased phason disorder. 

In the clean system, the local staggered magnetiza- 
tion has three different values, with the LSW results of 
m sl = 0.342,0.311 and 0.315 for z = 3,4 and 5, re- 
spectively. The averaged value determined by LSW is 
m s = 0.323, which compares well to the QMC result, 
m s = 0.3179(5) in the thermodynamic limit. Fig.[5]shows 
the averaged staggered magnetization plotted vs. the dis- 
order strength, as obtained by LSW and by QMC for a 
sample size of N — 1200 spins. The averaged value, 
rnj(A) increases with A for small A. This shows that 
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FIG. 4: (Color online) Plots of the integrated density of states 
N(E) as a function of energy (in units of J) for the pure 
system (blue) and averaged over the disorder for A = 0.025 
(red) for 1200 spins. 
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FIG. 5: Ground state staggered magnetization m s (A) plotted 
against A in LSW theory and QMC (sample size iV=1200). 



quantum fluctuations are reduced due to the presence of 
phason flip disorder, in keeping with the ground state 
energy dependence. This is a type of order-by-disorder 
phenomenon, and it can be contrasted with the effects of 
disorder in random bond and site- or bond-diluted an- 
tiferromagnets, where quantum fluctuations generically 
reduce the magnetic order. A possible exception was de- 
scribed in Ref. [22] for a doped quasi-one dimensional 
compound, although in that system the Neel tempera- 
ture is observed to fall with doping as well. To explain 
this manifestation of order-by-disorder a qualitative ar- 
gument goes as follows: it can be shown by considering 
a Heisenberg model on a cluster 15, [r| that the local 
staggered magnetization tends to be largest on sites of 
smallest coordination number. In effect, transverse spin 
fluctuations tend to be smaller when the number of near- 
est neighbors is smaller than the average value. Now, 
phasons in the pure SDL lead to the appearance of new 
z = 3 and z = 6 sites. The net change in m s is positive, 
due to the dominant contribution of small- z sites. This 
leads to the observed effect, namely an increase of the 
order parameter with phason disorder. It is interesting 
to compare the disordered staggered dice lattice with the 
results obtained in Ref. [23j for a disordered quasiperiodic 
antiferromagnet. In that case, the ground state energy 
is lowered - the quantum fluctuations get bigger as pha- 
son disorder is increased. In the quasicrystal, phason flip 
disorder causes magnon wavefunctions to become more 
extended and the density of magnons is increased, while 
in the staggered dice lattice, the situation is reversed. 

In conclusion, we have presented quantum effects aris- 
ing due to geometric disorder in a two dimensional an- 
tiferromagnet. The disorder does not introduce frus- 
tration, and it conserves the classical ground state en- 
ergy of the system. Contrarily to other disordered mod- 
els, phason disorder does not weaken antiferromagnetism 
which is, on the contrary, strengthened. The changes in 
ground state energy and the staggered moments indicate 
that quantum fluctuations decrease with increasing dis- 



order. This type of disorder could be realized in quasi- 
two-dimensional structures in which the nearest neighbor 
environments allow for independent local fluctuations be- 
tween two different conformations with little or no energy 
cost. We have discussed in this paper the hexagonal sys- 
tem, in which the two local bond configurations are iden- 
tical up to a reflection. The notion of phason disorder can 
be generalized to other systems such as the square lat- 
tice, where similar effects are expected to occur (work 
in progress). Realizations of this disordered model us- 
ing atoms in optical lattices, along the lines discussed in 
111 ], or by Josephson junction arrays can be envisaged in 
order to experimentally study the phenomena discussed 
here. 
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